Zeptosecond 7-ray pulses 
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High-order harmonic generation (HHG) in the relativistic regime is employed to obtain zeptosec- 
ond pulses of 7-rays. The harmonics are generated from atomic systems in counterpropagating 
strong attosecond laser pulse trains of linear polarization. In this setup recoUisions of the ionized 
electrons can be achieved in the highly relativistic regime via a reversal of the commonly deterio- 
rating drift and without instability of the electron dynamics such as in a standing laser wave. As 
a result, coherent attosecond 7-rays in the 10 MeV energy range as well as coherent zeptosecond 
7-ray pulses of MeV photon energy for time-resolved nuclear spectroscopy become feasible. 
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The time-resolved monitoring of fast-evolving pro- 
cesses with the pump-probe technique requires short laser 
pulses. The dynamics of chemical reactions has been 
probed by sub-picosecond laser pulses Strong in- 
frared laser pulses of femtosecond duration have been 
used to investigate the time evolution of vibrational wave 
packets in molecules 0- The new emerging techniques 
for generation of attosecond pulses of extreme ultravio- 
let radiation [sl and of recolliding electron wavepackets 
[3| have been exploited to look into the Auge-decay pro- 
cess 0, the nonlinear ionization dynamics molecu- 
lar dynamics Q and in theory to track the motion of 
an electronic wave packet in an atom Q. The required 
frequencies of the short pulses depend on the character- 
istic energies of the processes by means of which the fast 
dynamics is governed. While chemical reactions can be 
controlled with a few eV excitations driven by a weak in- 
frared laser field via one-photon processes, the molecular 
dynamics can be mediated with a few 10 eV ionization 
transitions by a strong infrared laser field via multipho- 
ton processes, and photon energies from a few 100 eV 
up to several keV are required to control the inner-shell 
electron dynamics. 

The time-resolved investigation of nuclear processes is 
a challenging problem How large photon energies 
and how short photon pulses are required to deal with 
this task? It is known [lO| that typical energies of nuclear 
single-particle transitions are of order of 1 — 10 MeV with 
typical decay lifetimes of the levels of around 10~^ — 10^^^ 
s. The energies of the collective nuclear excitations range 
from several dozens of keV up to 30 MeV; the electromag- 
netic giant dipole resonance is at about 15 — 22 MeV (the 
width of the resonance is of the order of 2 — 7 MeV) . This 
sets the scale for the required photon energies. The dis- 
integration time of the compound nuclei during nuclear 
reactions ranges from 10^^^ — 10^^® s. Some nuclear pro- 
cesses, such as the decay of excited levels, can be tracked 
with pulses longer than 100 fs for which the synchrotron 
radiation sources are well suited. However, there is a 
wealth of nuclear phenomena for which the investigation 



of the time resolved dynamics requires much shorter pho- 
ton pulses of up to a few zeptoseconds (zs) duration, such 
as e.g. resonance fluorescence [lH (1 fs timescale), reso- 
nance internal conversion [l^ (1 as timescale), compound 



nuclei evolution |13| or photodisintegration of nuclei 14 1 
(zs timescale). 

One of the successful ways towards coherent high fre- 
quency and ultrashort pulse production is connected with 
HHG faj. The highest photon energy achieved via HHG 
in gas jets today is at about 2.5 keV 15] and the shortest 
pulse length achieved in the same way is 130 as [l3|- To 
obtain higher photon energies, relativistic laser intensi- 
ties are required. However, the relativistic drift of the 
ionized electron suppresses the HHG yield in atomic sys- 
tems 17|. In the weakly relativistic regime there are 
methods to counteract the drift [l3| such that higher 
photon energies but not necessarily short coherent laser 
pulses are feasible. State-of-the-art proposals with zep- 
tosecond pulses reach keV photon HHG via overdense 
plasmas in ultrarelativistic laser pulses [Toj . 

In this letter we show that coherent 7-rays up to the 
10 MeV energy range as well as zeptosecond pulses of 
7-rays of MeV photon energy are feasible to allow for 
time-resolved nuclear spectroscopy. This is realized via 
relativistic HHG with highly charged ions in counterprop- 
agating strong attosecond pulse trains (APTs) of linear 
polarization (see Fig. [1]) . In this setup the common prob- 
lems with the relativistic drift of the ionized electrons 
at cutoff energy are circumvented: the electron is ion- 
ized by one attosecond pulse, driven by it up to the end 
of the pulse, then accelerated by the following counter- 
propagating pulse which reverts the drift and realizes the 
rescattering with the atomic core. 

We consider the HHG process of an atomic system 
driven by counterpropagating APTs which arise most fa- 
vorably from the same beam via a beam splitter. Our in- 
vestigation of HHG is based on the solution of the Klein- 
Gordon equation in the strong field approximation (SFA) 
[2^ . The APTs are linearly polarized with vector poten- 
tials Ai(a;) = A{tt + rii)ex and A2{x) — —A{ri2)ex, with 
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FIG. 1; With two counterpropagating APTs rather than a 
conventional sinusoidal pulse, it is possible to essentially re- 
vert the relativistic drift of the ionized electron and to enable 
electron rescattering in the highly relativistic regime (see the 
electron trajectory), x and z are the laser polarization and 
propagation directions of the APT, respectively. 
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FIG. 2: (a) Phase dependence of the electric field {E{rj)) of the 
counterpropagating APTs. The repetition angular frequency 
of the pulses in the train is tj = 0.05 a.u.. The propagation 
direction of the pulses is shown by the arrows; (b) Frequency 
spectrum (|c(fe)|) of one APT. The frequency components are 
phase locked. 



rji = t — z/c, 1^2 = t + z/c and the unit vector in the polar- 
ization direction Sx- The amphtude of HHG within the 
SFA in the single-active electron approximation is given 
by the following expression [2^ (atomic units are used 
throughout the paper): 



X Vh{x')G{x',x")Va{x"Mx")} 



(1) 



with harmonic order n, Vh{x) = 2Ah{x)/c ■ (p + 
Ai,(x)/c), VAix) = 2iV{x)/c^dt + V{xf/c^, atomic po- 
tential V{x), momentum operator of the electron p, vec- 
tor potential of the overall laser field in the radiation 
gauge Ai = A1-I-A2, matrix element of the vector poten- 
tial of the high harmonic field for an one-photon emission 
process Ah{x), bound state wave function <I>(a;) as an 
eigenstate of the energy operator in the radiation gauge 
24j . time-space coordinate x — {ct, x) — {ct, x, y, z), and 
speed of light c. In Eq.ijT]) G{x',x") represents the Green 
function of the Klein-Gordon equation for the electron in 
the field of both APTs. There are two different scenarios 
for the electron dynamics after ionization. In one sce- 
nario the electron successively moves in different coun- 



terpropagating pulses. In the other scenario, the two 
counterpropagating pulses act on the ionized electron si- 
multaneously during its excursion time. The recoUisions 
of the ionized electrons with non-vanishing probability 
and HHG are connected only with the first scenario of 
the interaction in the highly relativistic regime because 
of chaotic electron dynamics and a negligibly small range 
of electron phases for rescattering in the second scenario 
2l|. We thus restricts ourselves to situations where the 
first scenario applies. Without loss of generality the pro- 
cess is initiated by the action of the first laser field fol- 
lowed by the second one. Then, the Green function of 
the Klein-Gordon equation for the electron in the field of 
both APTs can be approximated via the Klein-Gordon 
Volkov Green function for a single laser field GY{x',x") 
(i G {1, 2} refering to the respective APT) in the follow- 
ing way: 



G^{x,x') = i 



d^^^G^{x,x^)d,t. 



Gr(x^,x') (2) 



with time such that the first laser pulse has already 
left the wave packet of the active electron and the second 
laser pulse has not acted yet on the electron at space 
coordinate x^. As in the parameter range of interest 
here Klu <C Ip with the largest frequency component of 
the attosecond pulse K, repetition angular frequency in 
the APT to and ionization potential Ip, the integral in 
Eq.ijT]) can be calculated via the saddle-point method. 

In Fig. [3] the HHG yield is shown in the setup of coun- 
terpropagating APTs for a pulse shape as in Fig. [2l The 
yield is compared with the spectrum calculated in the 
dipole approximation (DA) for the setup of two coun- 
terpropagating APTs as well as with the result of the 
specially tailored APT as in [2^ with pulses of rectan- 
gular shape. The DA result can be regarded as a refer- 
ence for the ionization-rescattering process without drift, 
since in this description the relativistic drift is omitted. 
In the considered strongly relativistic regime with a laser 
field amplitude of Eq = 88 a.u. (peak intensity 6 x 10^° 
W/cm^), the counterpropagating APT setup produces 
harmonics with high efficiency, matching the DA result, 
in the cutoff region of 1 MeV energy. It shows also a 
strong improvement with respect to the proposal which 
employs a single tailored APT gaining up to three 
orders of magnitude in the HHG rate. Moreover, the con- 
ditions for the pulse features in the present setup are less 
demanding as those for the tailored pulses in [20| . In the 
setup here each pulse in the train is of about 300 as du- 
ration and contains less than 20 harmonics, while for the 
single tailored APT up to 100 harmonics are required. 
Note that the severe requirements for the special pulse 
tailoring is absolutely necessary in the case of a single 
APT as the negligible HHG yield in the pulse without 
tailoring in Fig. shows. 

The dependence of the HHG spectra on the pulse sep- 
aration is presented in Fig. [3][b). The HHG spectra are 
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FIG. 3: Harmonic emission rate in APT propagation direc- 
tions as function of the harmonic energy {dwn/dQ. is the emis- 
sion probability of the n"* harmonic per unit time and per 
unit solid angle fi). The laser field amplitude is Eq — 88 a.u., 
Ip = 63 a.u. (Mgi°+): 

(a) (solid black) two counterpropagating APTs via the Klein- 
Gordon equation, (solid gray) two counterpropagating APTs 
within the DA, (dotted) tailored APT as in |20| via the Klein- 
Gordon equation, (dashed) two copropagating APTs with the 
pulses as in Fig. [2] via the Klein-Gordon equation (see very 
bottom left), uj = 0.05 a.u.. The window shows enlarged the 
cutoff energy region, (b) The dependence of HHG yield on 
the time delay between the pulses; the delays are indicated 
for each spectrum in units of fs. 



characterized by a bent plateau that ends in a sharp cut- 
off. The bending increases for shorter time delays be- 
tween the pulses, whereas the oscillation pattern due to 
the interference of two possible quasiclassical trajecto- 
ries decreases. Both features indicate stronger relativis- 
tic signature of the process for shorter time delays [1^ . 
The longer the time interval between the two pulses, the 
more excursion space is given to the electron and the rec- 
oUiding electron need to ionize at higher field strengths 
resulting in higher ionization and thus higher HHG rate. 
Additional quantum spreading of the ionized electron can 
reverse this effect for long time delays as visible from the 
spectrum. Further, the cutoff energy is slowly increas- 
ing with a reduction of the time delay. This is because 
the ionization of the recoUiding electron tends to occur 
at a later moment in the first pulse if the time delay is 
shorter. As a result, due to the particular shape of the 
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FIG. 4: (Color online) Harmonic emission rate in APT propa- 
gation directions as function of the harmonic energy for coun- 
terpropagating APTs with various time delays compared with 
the spectrum in the DA (pulse delay 10 fs) (gray). The pulse 
delays for the corresponding spectrum are indicated in units 
of fs; further, Eo = 325 a.u. and Ip = 150 a.u. (Ari^+). 



pulses in Fig. [21 the energy of the electron after the first 
pulse is larger and the time left for the final acceleration 
in the second pulse is higher, too. Both dependences, the 
increase of the ionization rate as well as the decrease of 
the cutoff energy become weaker for longer time delays. 

The counterpropagating APTs setup enables the ex- 
tension of efficient HHG into the highly relativistic 
regime. Thus, in Fig. [4] we show HHG spectra employ- 
ing APTs with Eo = 325 a.u. (peak intensity 8 x 10^^ 
W/cm^) and various time delays between the APTs. The 
harmonic emission rate in the cutoff region is as intense 
as the one in the DA, i.e. the recoUision can be realized 
even in the highly relativistic regime. Then, HHG cutoff 
energies of about 10 MeV with a significant emission rate 
are possible. The rate decrease compared with the previ- 
ous case of Fig. |3]is only due to quantum spreading which 
is enhanced because of the stronger laser electric field and 
the longer time delay between the pulses. The latter is 
necessary for reverting the relativistic drift. Further, the 
reduction of the rate at cutoff energy, when decreasing 
the time delay from 10 fs, is larger in this regime than 
for weaker laser fields. The cutoff energy dependence on 
the time delay is inverted compared to the weak field case 
due to the stronger relativistic signature of the process. 
An optimal time delay is at about 8 fs. 

In the relativistic regime, a broad spectral content of 
the emitted harmonic radiation as well as the fact that 
mainly one trajectory contribute to HHG, are beneficial 
for ultrashort pulse generation. An analysis of the har- 
monic phases shows that the variation of the rescattering 
phase by 0.35 yields a variation of the rescattering energy 
of 1000 keV. Then the spectral window to generate the 
shortest possible pulse is at about a few keV. Choos- 
ing a 4 keV window near the spectral region of 1 MeV 
where the HHG yield is largest, a zeptosecond pulse with 
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FIG. 5: Temporal shape of the generated pulse with a Gaus- 
sian window at 805 keV with a FWHM of 4 keV. The FWHM 
of the pulse is approximately 900 zs. Eq — 88 au., Ip = 63 
a.u. and a time delay between the pulses of 0.9 fs. 



FWHM of 900 zs can be produced (see Fig. [5]). The 
zeptosecond pulse in the considered case is rather weak: 
= 10^ coherent photons per pulse are feasible when an 
interaction length of 1 cm, a cross-section of 1 mm^ and 
an ion density of 10^^ cm~'^ (e.g. in the "bubble" regime 
of laser-plasma interaction [25|) are employed. However, 
due to the shortness of the pulse, the peak flux density of 
the photons is rather large: 10^^ photons- s~^mm~^ |2q . 
On roughly estimating the number of interaction events 
N — aNyPnLint between zeptosecond pulse and nuclei, 
with cross-section a of photonucleus collective interac- 
tion, nuclei density p„, and interaction length Lint, we 
obtain N ~ 10^ interaction events per zeptosecond pulse 
with a ~ 10-26 cm2, p„ ~ 10^3 cm'^, and Lint ^ 10^ 
cm. The latter indicates at least theoretical feasibility 
of nuclear spectroscopy with these pulses. We note that 
the intensity of the zeptosecond pulse could be increased, 
e.g. by an additional tailoring of the pulses in the APTs. 

In conclusion, we have shown that HHG in the rel- 
ativistic regime with an atomic system can be realized 
employing counterpropagating strong APTs. Coherent 
7-rays in the 10 MeV energy range as well as coherent 
zeptosecond pulses of MeV photon energy can be gener- 
ated this way. The main difficulty in experimental real- 
ization of this setup is related to the availability of strong 
APTs. Those pulses combine two extreme properties: 
short pulse duration and high intensity. However, recent 
advances in the field of strong laser field interaction with 
overdense plasmas show that strong attosecond pulses are 
expected to be created in the near future [ll,[23. 
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